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Abstract. For monotone linear differential systems with periodic coefficients, the (first) Floquet eigen- 
value measures the growth rate of the system. We define an appropriate arithmetico-geometric 
time average of the coefficients for which we can prove that the Perron eigenvalue is smaller 
than the Floquet eigenvalue. We apply this method to Partial Differential Equations, and we 
use it for an age-structured systems of equations for the cell cycle. This opposition between 
Floquet and Perron eigenvalues models the loss of circadian rhythms by cancer cells. 

Une inegalite pour les valeurs propres de Floquet et de Perron de systemes 
differentiels monotones et d 'equations structurees en age 

Resume. La (premiere ) valeur propre de Floquet decrit le taux de croissance des systemes differentiels 
lineaires monotones a coefficients periodiques. Nous definissons une moyenne arithmetico- 
geometrique en temps des coefficients, qui nous permet de demontrer que la valeur propre 
de Perron pour le systeme ainsi moyenne est plus petite que celle de Floquet. La methode 
s 'applique aux Equations aux Derivees Partielles et nous I ' utilisons pour un systeme d' equa- 
tions structurees en age qui decrit le cycle cellulaire. Cette opposition entre valeurs propres 
de Floquet et de Perron modelise la perte de controle circadien pour le cycle cellulaire des 
cellules cancereuses. 



Version fmncaise abregee 

Les systemes biologiques sont souvent soumis a des controles periodiques. Un exemple en est fourni 
par le rythme circadien (journalier) qui trouve son origine au niveau des noyaux suprachiasmatiques de 
F hypothalamus dans un reseau de regulation genique a present bien etudie J9][l]. Une question medicale 
reliee est de comprendre l'influence de ce rythme sur la croissance des populations de cellules tumorales, en 
particulier dans une perspective therapeutique. Un cadre mathematique naturel pour etudier cette question 
est fourni par la theorie des equations physiologiquement structurees. 



The second author was partially supported by the joint RFBR-CNRS grant number 05-01-02807 



J. Clairambault, S. Gaubert, B. Perthame 



Dans cette note, nous proposons de comparer la (premiere) valeur propre de Floquet, qui decrit la crois- 
sance d'un systeme dont les parametres sont soumis a un controle periodique, et la valeur propre de Per- 
ron, qui decrit la croissance du meme systeme, mais a coefficients moyennes. Pour parvenir a cette com- 
paraison, nous introduisons une moyenne arithmetico-geometrique appropriee, dans trois cas : Equations 
Differentielles Ordinaires, systeme dynamique discret, Equations aux Derivees Partielles (EDP) structurees 
en age. 

Le premier cas concerne un systeme differentiel lineaire monotone periodique. Soit t i— > A[t) une ap- 
plication T-periodique, a valeurs dans K dxd , integrable sur [0, T], et considerons l'equation differentielle 
X(t) = A(t)X(t), oil X est absolument continue et X(0) est present. Nous supposerons que pour i ^ j, 
on a Aij (t) > pour presque tout t, ce qui garantit que le flot en temps positif associe a cette equation 
differentielle preserve l'ordre partiel usuel de M. d . Les proprietes spectrales de cette equation differentielle 
relevent alors de la theorie de Perron-Frobenius. En particulier, la premiere valeur propre de Floquet, A per , 
est le plus grand reel tel qu'il existe une fonction absolument continue X, T-periodique, a valeurs dans 
(ou R + designe Fensemble des reels positifs ou nuls), et non identiquement nulle, telle que l'equation 
differentielle (fTJ soit satisfaite. Rappelons que la moyenne arithmetique d'une fonction T-periodique u(t) 
est donnee par (fJJ, et definissons la matrice A par ©. Cette matrice a coefficients hors diagonaux positifs 
ou nuls est dotee d'une valeur propre de Perron classique A s , qui est le plus grand reel tel qu'il existe un 
vecteur non-nul U £ tel que \ S U = AU. Le resultat suivant peut etre vu comme une generalisation de 
celui de Q. 

THEOREME 1. - On a toujours A per > A s . 

Une inegalite analogue est aussi obtenue pour les deux autres cas considered, d'un systeme discret en 
temps et du systeme d'EDP structurees en age © qui decrit le cycle de division cellulaire. 



I Introduction 

Biological systems are often subject to periodic controls. This occurs for instance with circadian rhythms, 
the origin of which are found in the suprachiasmatic nuclei of the hypothalamus, in a now well established 
gene regulatory network J9j [T] . A related medical question is to understand the interactions between the 
cell cycle and this circadian rhythm, which is expressed in every nucleated cell, with coordinating inputs 
from the hypothalamus. How can these rhythms induce differentiated growth between healthy and tumoral 
cells? A molecular mechanism has been evidenced [9], but most importantly in laboratory experimental 
settings, tumour growth has been shown to be favoured by disruptions of the normal circadian rhythm H. 
From a mathematical point of view, cell population growth is well described by physiologically structured 
equations, see [7, 2], and the first eigenvalue of the underlying differential operator is the natural quantity 
that accounts for the growth of the system. Our purpose is to compare the first eigenvalues in the case 
of constant and periodic coefficients. It is simpler to consider in the first place differential systems; with 
periodic and nonnegative coefficients this eigenvalue is nothing but the (first) Floquet eigenvalue; with 
constant coefficients it refers to the usual Perron eigenvalue. Surprisingly, it is possible to prove in great 
generality that the Floquet eigenvalue is larger than the Perron eigenvalue with an appropriate arithmetico- 
geometric average of the coefficients. Precise statements, and proofs, are given in the first subsection in the 
case of a differential system. The result also extends to discrete time systems (section ITITb and to Partial 
Differential Equations and we give, in a fourth section, an application to age- structured systems. Finally in 
a fifth section, we briefly comment on the relevance of these models to physiological systems. 

II Differential systems 

Let t I— > A(t) be a T-periodic map with values in M. dxd , integrable on [0, T], and let us consider the 
differential equation X{t) = A(t)X(t), with a prescribed initial condition X(0) £ R d (when dealing with 
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such differential equations, we will always require X to be absolutely continuous, and we understand that 
the equality holds for almost all t). We assume that for i ^ j, we have Aij(t) > for almost every t, 
so that the flow in positive time of this differential equation preserves the standard partial ordering of M d . 
Hence, the spectral properties of this differential equation belong to Perron-Frobenius theory. In particular, 
the (first) Floquet eigenvalue, A per , can be introduced by means of the following positive Floquet problem: 
there exists a T-periodic function X, with values in El (E + denotes the set of nonnegative real numbers), 
non identically zero, such that 

X(t) = A{t)X{t) - X per X{t) , t e R . (1) 

If A satisfies some irreducibility properties, for instance if Ay (t) > for almost every t and for i ^ j, A per 
is uniquely defined by the previous property, and X is unique up to a multiplicative constant. However, the 
results of this note are also valid in the reducible case, in which the Floquet eigenvalue can be defined as 
the maximal real number A per such that there exists a function X with the above properties. 
We now introduce the arithmetic mean of a T-periodic function u(t) as 

1 " T 



(u) a = -J u(s)ds (2) 

and we define the constant coefficient matrix A with entries: 

An = (A u ) a , l<i<d, A i:j =exp((log(Ay))J i^j,l<i,j<d. (3) 

Observe that (log(Ay)} S K U {— oo} is well defined for i ^ j, because the positive part of log Ay(t) is 
integrable as soon as A^ (t) is integrable. Since the off diagonal coefficients of the matrix A are nonnega- 
tive, we can apply the Perron-Frobenius theory and consider its first (Perron) eigenvalue, X s , which is the 
maximal real number such that there is a non-zero vector U E such that X S U = All. As discussed in 
SectionM the following result generalises the one of 0. 

Theorem 2.1. - We have always X per > X s . 

Proof. - We first prove the inequality when Aij(t) > 0, for almost all t and for i ^ j. Then, the 
function X in (Q]i is such that Xi(t) > for all t and for all i. We set, 

Xi (t) := log X 2 (t) and for i ^ j, l l} (t) := log A l} (t) . 

From the differential system (Q~|), we obtain 

±i(t) = J2xrHt)A ij (t)X j (t) - A per 
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eXp(-Xi{t) + Uj{t) + Xj[t)) + A u (t) - A per 



Taking first the arithmetic mean on [0, T] componentwise, and then using Jensen's inequality, it comes for 
all 1< i < d, 



0= (J2exp(-Xi(t) + iijit) + Xj (t)) ) +(A u (t)) a -X 



per i 



> J2 e M~ (Xi(t)) a + (l l3 (t))a + <*i(*)> tt ) + (Mtya - V 
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Setting Xi := exp((xi(t)) ), with the definition ([3]) of A, this inequality also reads 

> 1 A i] Xj - Aper , 



and thus, multiplying by Xi, 



AX < X pa X 



Using the Collatz-Wielandt characterisation of the Perron eigenvalue of A, X s = min{r; 3Y £ intR+,Ay < 
rY}, we deduce that X s < A per . The general case is obtained by considering the matrix with entries 
Aij (t) + e, with e > 0, and by applying a continuity argument, the details of which are left to the reader. □ 

III Discrete systems: an inequality for the Perron eigenvalue of a geometric mean 

The same proof allows us to treat discrete systems. For k E N, let A(k) = A(k+p) beadxd, p-periodic 
matrix with nonnegative coefficients. We define A per to be the maximal real nonnegative number such that 
there exists a non-zero p-periodic solution X(k) with values in to 

X pei X{k+l) = A{k)X(k) , (4) 

so that Aper is the Perron eigenvalue of the product A(p — 1) ■ • • j4(0). We now set (u) a :— p _1 (u(0) + 
• • • + u(p — 1)) for all p-periodic functions u, and we define the constant coefficient matrix A, with entries: 

A i:j = cxp ((\og(A tJ )) a ) , 1 < i,j < d (5) 

(unlike in the continuous time case, the nature of the mean is the same for diagonal and off diagonal entries). 
We denote by A s the Perron eigenvalue of the nonnegative matrix A. 
Theorem 3.1. - We have again Ap e r ^ A s . 

Proof. - Since the Perron eigenvalue is a continuous function of the entries of a matrix, it suffices to 
show the inequality when the entries Ay(fc) are all positive. Then, the vectors X(k) above have positive 
entries. We set 

Xi(k) := logXi(k), iij(k) := log Aij (k), ^ := log(A pei ) . 
We take componentwise the logarithm in X pel X(k + 1) = A(k)X(k) and arrive at 

Xi(k + 1) -Xi(k) =log |^exp ( - Xi(k) +£ij{k) + Xj(k)) j - p . 
Taking the arithmetic mean for k = 0, . . . , p — 1, it comes 



= / log ^exp ( - Xi(k) + £ij(k) + Xj(k)) 



A 4 • 



Next, because the function f{y\, . . . , yd) = log(J^ - exp(?/j)) is convex, we apply Jensen's inequality and 
obtain 

0>log (^e^{-{ Xi {k)) a + {iij{k)) a + {xj{k)) a ^ -fi , 
and, exponentiating, we get 

x pei x > Ax , 

with Xi := exp((xi(k)) a ). Using again the Collatz-Wielandt characterisation of the Perron eigenvalue of 
A, we deduce that A per > A s . □ 
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IV An age-structured system for the cell division cycle 

General references and experimental validations on the topic of structured population dynamics and cell 
cycle can be found in 02] |7). For a recent mathematical approach based on entropy properties, we refer 
to (U [lO). Here and following earlier work [3|, we model our population of cells by a Partial Differential 
Equation for the density rii(t, x) > of cells with age x in the phase i = 1, . . . , at time t, 

{ gj«i(*> x) + ifcni(t, x) + [di(t, x) +Ki-> i+1 (t, x)}rii(t, x) = , 



rii(t, x = 0) 



JQ_i^i(t,x') m-i{t,x') dx', 2 < i < I, 



x'>Q 



ni (t, x = 0) = 2 / Ki—,x(t> %') nj(t, x) dx' 



(6) 



Here and below we identify I + 1 to 1. We have denoted by di(t, x) > the apoptosis rate, by K^i + i the 
transition rates from one phase to the next, and the last one (i = I) is mitosis where the two cells separate. 
These coefficients can be constant in time (no circadian control) or time T-periodic in order to take into 
account the circadian rhythm. Our assumptions are 



K 



(t,x) > 0, di{t,x) > 



are bounded, 



(7) 



and, setting 



min Ki^ i+ x(t,x) := ki^ i+ x(x), max [di + K^+i] := ^i(x), M, ; (x) = / Hi{y)dy, 



n 



(j/)e' 



-Mi(y) 



dy > 1/2 



(8) 



With these assumptions and following [8|, one can again introduce the growth rate (Floquet eigenvalue) of 
the system: A pe r E M such that there is a unique T-periodic positive solution to the system: 

' £ t Ni(t,x) + &Ni(t,x) + [di{t,x) + A per + Ki^ i+ x(t,x)}Ni(t,x) = 0, 



Ni(t,x = 0) 



a/) iVi-i^x') dx' , 2 < i < I , 



z'>0 



(9) 



Nx(t,x = 0) = 2 [ Ki^x{t,x')Ni{t,x')dx' / JV 4 (t,x)dx = l 

As in formula ((3), we can define the averages 

| (di(x)) a = tj; [ di(t,x)dt, (Ki^ i+ x(t,x)) a = i / Ki^ l+1 (t,x)dt , 



(Jfi_>i+i(t,x)) = exp 



log (if^j+i^x))^ 



(10) 



These averages define the Perron eigenvalue A s G R, which is such that there is a unique positive solution 
to the system 

< &Ni(x) + [(d i (x)) a + \ s + (K i ^ i+1 (t,x)) a ]N i = , 



Ni(x = 0) 



(K i - 1 ^ i (t,x')) a N i - 1 (x>)dx', t^l , 



x'>0 



(ID 



A^x = 0) = 2 / (X / ^ 1 (i,x'))„/Y/(x , )rfx' 

Jx'>0 
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We have the following analogue of Theorem 12. II 

THEOREM 4.1. - Under assumptions 0-(®, we still have A per > X s . 
Observe that the arithmetic mean of the coefficients K^i+i is taken in the PDEs, whereas their geometric 
mean is taken in the integral equations. Hence, an artificial loss rate of cells (l"Q_>j_|_i) a — (Ki^i + i) g from 
phase i to phase i + 1 arises in the averaged model. 

Proof. - The proof differs slightly from that of Theorem 12. 11 because, working with x 6 R + the 
corresponding quantities cannot always be normalized as measures. Therefore, following [3|, we define 
qi(x) = (log Ni(t, x) — logiVj(a;)) and we have, up to the insertion of a factor 2 when i = 1, 

f n\ A f^MA A „ [/v , vv iVi-i(s) Ki-i^x) , \ 

''^^rw'/rri <*-x-*(»)>. ^ -——-^^ . 

We can now define the probability measures 

d ^(x)=(K i _ 1 ^ i (x)) g ^tM fori ,4 1, ^(^=2 (^(^^H . 
Using Jensen's inequality, we obtain 

= > ( / log -?y — t^7t? rrH^MiW > 

\J Ni-iO) (Ki-x^i(x)) g I a 

JVi(M) \ / / ^_i^(t,a;) , A 



JVi-i(a:) 7„ V (AVi^W) 

qi-i(x)d(i,i(x) 



since by definition of (l£i_i_^(a;)) we have / / log( ' 1 ~ > '( ) \ d/ij (x) = 0. Because qt satisfies 



— <7i + Aper - A s = , 



8_ 

dx 

we arrive at 

qi(x = 0) > J qi-i(x) dfii(x) = J [qi-i(0) + (\ s - \ pel )x] dpi(x) 
Therefore, summing over i from 1 to /, we have obtained the result since 
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> (A s - Aper) ^2 X dfli(x) 
»=1 ^ 



□ 



V Physiological relevance of the model to the question of circadian control on tumour growth 

In 0, with different averaging of transition rates Ki—*i+\, we had found no natural hierarchy between 
Aper and X s , whereas periodic control exerted on the sole apoptosis rates di produced a Floquet eigenvalue 
which was shown to be always higher than the corresponding Perron eigenvalue obtained by arithmetic 
averaging of the rfj. The present note gives natural hypotheses that are needed to compare the eigenval- 
ues and conclude that Ap er > A s when the coefficients Ki-^i+i also vary. From these results a question 
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arises: if A per > X s , then how come that experimental results [4] show faster tumour growth when the 
normal circadian rhythm is disrupted by irregular light inputs? Experimental tumour growth curves with 
and without jet-lag-like circadian clock disruption J4[ essentially differ at the beginning of proliferation, 
when it shows almost pure exponential behaviour, so that possible discrepancies between theoretical and 
experimental results are most likely not due to not taking into account nonlinear feedback. Clearly, setting 
at a constant value (instead of periodic) the coefficients di or Ki^i + i is not sufficient to take into account 
the complexity involved in the disruption of circadian control. To this purpose, it may be necessary to 
represent more complex control mechanisms involving, e. g., the inhibition of cyclin dependent kinases by 
clock-controlled genes such as Weel together with elaborate models of the molecular circadian clock 
||6l and its disruptions. 
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